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Abstract 

We present an analysis of the continuum extrapolation of fs in the 
static approximation from lattice data. The method described here aims 
to uncover the systematic effects which enter in this extrapolation and 
has not been described before. Our conclusions are that we see statistical 
evidence for scaling of /^*"* for inverse lattice spacings ^ 2 GeV but 
not for ^ 2 GeV. We observe a lack of asymptotic scaling for a variety 
of quantities, including /^*"*, at all energy scales considered. This can 
be associated with finite lattice spacing systematics. Once these effects 
are taken into account, we obtain a value of 230(35) MeV for Z^*"^* in 
the continuum where the error represents uncertainties due to both the 
statistics and the continuum extrapolation. In this method there is no 
error due to uncertainties in the renormalization constant connecting the 
lattice and continuum effective theories. 



1 Introduction 



The decay constant of the S— meson, fs, defined through the matrix element of 
the axial current, A^, 

<0\A,\Bip) >=tfBP>. (1) 

(where B{p) is a 5— meson of momentum p), is an essential ingredient in many 
calculations in the Standard Model of Particle Physics. For example, it enters in 
the theoretical determinations of (i) the mass splitting in the B — B system, which 
is proportional to the square of fs and, in principle, (ii) the Cabibbo-Kobayashi- 
Maskawa matrix element, Vub-, and hence, through the unitary triangle to the 
element Vtd (see eg. IH). Clearly therefore, for a more complete knowledge of 
some of the fundamental parameters of the Standard Model, theorists need to 
provide a prediction of fs (see also 0). 

The two main methods of determining Jb from theory are QCD Sum Rules 
(see eg. 0) and Lattice Gauge Theory (see eg. @, 0). The QCD Sum Rules has 
an important advantage compared to the lattice in that it is an analytic, rather 
than a numerical method. However it involves approximations to the full theory 
that cannot be systematically improved, and it has a perturbative component 
calculated at a momentum scale ~ 1 GeV where the validity of perturbation 
theory can be questioned. 

In contrast to QCD Sum Rules, the lattice is a fully non-perturbative ap- 
proach and furthermore uses approximations that can be systematically im- 
proved. For each approximation there is an associated tunable "lattice parameter" 
which can eventually be adjusted towards its physical value thereby removing the 
approximation. For example, lattice calculations are performed on a finite phys- 
ical volume which can, in principle, be increased towards the physical value (ie. 
infinity) such that finite volume effects can be neglected to within any desired 
accuracy. In general the size of the systematic errors introduced due to these 
lattice approximations has been quantified and often does not present a problem 
given the size of the statistical errors. This is particularly true in the following 
analysis. 

A lattice approximation that warrants further discussion is the "Quenched 
Approximation" where the effects of sea quark loops are suppressed. All the 
lattice results discussed here were obtained using the quenched approximation. 
There is, however, strong circumstantial evidence that its effects on hadronic 
masses are small and its effects on hadronic matrix elements are of the level 
of the statistical errors (10% - 20%) There will be a further discussion on 
quenching in the final section. 

The reason that these lattice parameters are not set to their physical values in 
the first place is due to the limitations of present computing power. In this sense, 
with the inevitable development of computing resources, lattice QCD results will 
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be able to provide predictions with ever increasing precision. 

Turning to the determination of Jb, there is one immediate, technical problem 
in simulating a 6— quark on the lattice: its mass is larger than the ultraviolet cut- 
off provided by the finite spatial resolution of the lattice, ie. > a~^. Here a~^, 
the inverse lattice spacing, is typically of order 1-3 GeV in present simulations. 
There are two approaches to overcome this problem. One is to study quarks which 
are fairly heavy, but whose masses are still less than and then to extrapolate 
to rrif,, and the second is to use the "static approximation" which is the leading 
term in the heavy quark expansion [} This paper deals with lattice results 
using this second approach. To obtain the real value of fs (ie. to all orders in the 
heavy quark expansion), one must interpolate results from the two approaches 

i- 

Historically, lattice calculations of Z^*'^* (ie. J'b in the static approximation), 
have proved very difficult due to the low signal to noise ratio of the hadronic cor- 



relators required for the calculation 0. The first successful measurements [|T0|, [11 
produced surprisingly large values, f^""^ ~ 300 MeV. Since then, many calcula- 
tions, at the same value of a, have generally confirmed these early measurements 
[0-[|l9]. However when calculations at different values of a were performed there 
appeared, in some analyses, to be a trend towards smaller values of /^*"* as a — > 0, 
ie. as the "lattice parameter", a, is tuned towards its physical value p^ . 
The requirement for lattice results at a non-zero value of a to be physically rele- 
vant is that physical quantities remain a constant in the continuum limit, a ^ 0, 
ie. that they scale Q This apparent downward trend of /^*"*, if actually present, 
is therefore a violation of scaling and is to be contrasted with the apparent scaling 
behaviour observed in other physical quantities [^. It must be stressed that, 
statements regarding scaling or the lack of it can only be made within a certain 
statistical uncertainty, and for quantities as difficult to determine as /^*"*, these 
errors are typically large. At some point scaling violations may disappear below 
the level of statistics, and hence become unimportant. 

The issue of scaling is normally studied by plotting the ratios of physical 
quantities as a function of a, as was the case in the above-mentioned analyses of 
j^at which an apparent violation of scaling was observed. This paper outlines 
an alternative method of studying scaling behaviour in which the functional form 
of the dimensionless lattice quantity corresponding to /^*°*, (ie. Z^,, defined later 
in eq.(p!OD), is fitted as a function of the bare lattice coupling. We show that this 
analysis applied to the currently published data on /^*"* shows that its scaling 
violations for a"^ ^ 2 GeV are small, ie. that they have fallen to below the 
statistical errors. However, for ^ 2 GeV we do find scaling violations. This 



^The static approximation will be discussed in more detail in the next section. 
^This nomenclature is in analogy with critical phenomena where all physical quantities are 
proportional to the correlation length raised to a scaling index. 
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suggests that the conclusions of the studies |]5|, 0, |T3[ are due to the inclusion 
of data too far from the continuum, ie. with a^^ < 2 GeV. The implication of 
this result is that, at present, the best estimate for the continuum value of f§°'*' 
is given by those simulations at pa 2 — 3 GeV. 

In the following section we outline the static theory as applied on the lattice. 
Those readers familar with the static theory on the lattice may wish to skip this 
section. In sec. 3 we study the scaling of /^*"* and other physical quantities, and 
in sec. 4 we discuss our results. We also investigate finite lattice spacing, or 0(a), 
effects and the approach to asymptotic scaling. Once these issues are studied and 
understood we obtain a continuum value for /^*°* of 230(35) MeV. 



2 Static theory on the lattice 

This section begins with a brief description of the lattice method of calculating 
physical quantities. For a full discussion see, for example, the review articles 

Lattice QCD typically studies Euclidean correlation functions, C(t), of hadronic 
operators O^, 

C{t) =< 0^{t) 0^(0) > . (2) 
For purposes of explanation, we will define 

0^(t) = Y.M^.t) = E^i(5.^)7o75V'2(5,t), (3) 

X X 

where the 70,5 are Dirac spinor matrices, and the subscripts 1, 2 are flavour indices. 
At large times t, C{t) has the following behaviour: 

_^ I <0|O^|P> P Mj, t 

where P is the lowest state with the quantum numbers of O^, in our case a pseudo- 
scalar meson. Once the exponential behaviour in eq.(0) has been established, a 
simple exponential fit provides lattice predictions of the decay constant, 

/J; =< 0|O^|P > /M^, (5) 

(cf. eq.(|I|)) and the mass, Mp. We will refer to these lattice predictions, fp 
and Mp, generically as f2f (whether they be hadron masses or matrix elements). 
These are not physical quantities because (i) all the fields in the lattice action are 
made dimensionless using the lattice spacing a, and therefore all the predictions, 
f2f , from the lattice are also dimensionless, and, (ii) in the case of matrix elements, 
the operator O^, is not correctly normalized. 
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To determine the scale a we choose a lattice prediction of one physical quan- 
tity and compare it to the experimental (dimensionful) value, The lattice 
operator is correctly normalized to its continuum counterpart O typically 
through a multiplicative renormalization factor Z^*^", ie. O = Z^^" (see eg. 



2^). Combining both steps we have: 

1 ^7 



(6) 



(In this discussion we are assuming that f2j has energy dimension unity.) For 
non- matrix element quantities such as hadron masses M^, Z^^"' = 1. 

In lattice QCD, we have dimensional transmutation in action: we first set 
the bare coupling in the lattice action, g'^, at the start of the calculation, and 
then determine the corresponding ultraviolet cutoff, a~^, through eq.(|^). (We 
could instead proceed in the reverse direction, by first setting a~^, but this would 
require expensive simulations at many trial values of before we settled at the 
correct value of g"^ corresponding to the chosen value of a~^.) 

Finally then we have the lattice prediction of a dimensionful physical quantity, 



= Zf^'^aT' (7) 



The correlation function Cif) is numerically calculated using the Wick con- 
tracted form of eq.(ll), 

C{t) =< rr{Gi(t, 0)7075^2(0, t)7o75} > • (8) 

The quark propagators Gi^2 are defined from the lattice version of the Dirac 
equation. 

In the static version of the theory Gi is defined using the solution of the Dirac 
equation in the limit of an infinitely heavy quark (see eq.(2.37) in ^). In this 
case we are simulating a pseudo-scalar meson (which we'll denote B) made up of 
a light quark and an infinitely heavy quark. We have from eqs.(|l],§). 



where Eb is the binding energy^ of the light quark in the B-meson, and Mb is 
the experimental value of the 5— meson mass. We define 

zl^(n^^f=^J£-l^. (10) 



^ We have factored out the exponential time dependence of the (infinitely-heavy) quark in 
the static propagator. 
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This is the lattice quantity corresponding to f^'^^ (see sec.l). Note that in |]T2 
the symbol Jb is used for Z^. Correctly normalizing the operator, and putting 
in the appropriate dimensions (cf. eq.(0)) we have, 

The issue of scaling is then simply a question of whether the dependences of 
the Zj^^^, Zl and a~^/^ cancel in eq. (|Tl|) . 




Using eq.(|^) we have finally 



IB 



Mb (Zf^" nf')3/2- 



(12) 



3 Continuum scaling of f^^"^ 

The usual method of determining whether the scaling^ region of a lattice simula- 
tion has been reached is to study a dimensionless ratio, of physical quantities 
fli, to see if this ratio is a constant in a, or equivalently in g'^. Using eq.(0) 
the Qi in the ratio can be replaced with the dimensionless quantity . 
Dimensionless ratios of lattice quantities are generally used wherever possible be- 
cause statistical fluctuations, and indeed, some systematic effects tend to cancel. 
For this reason ratios are often better determined than the absolute quantities 
themselves. 

Using this "ratio" method, we show in fig.(|l|) the plot of /^*"* versus a^, ie. 
the lattice spacing determined from the string tension a. The data plotted is a 
collection of published data on /^*°* and is reported in table with references 
cited in column 1. (Note that P = G/g"^, where g is the bare lattice coupling 
which appears in the lattice action.) For completeness, in table |I], we also list 
the values of and used and the f^""^ obtained by each group. In the 
figure we use eq.([T^) to determine the /^*"* values given the Zl values in table |I[ 



The values we have used are the boosted, tadpole improved values 
which means, for instance, at /3 = 6.0 we have -Z^g" = 0.70 for the Wilson action, 
and 0.79 for the clover action. 

It is difficult to interpret this plot, and difficult to gauge the likelyhood of a 
decrease or increase of /^*°* as a — 0. Of course, one can simply use a linear or 
quadratic fit 0, ^, 0, |l3l, and fit the data in any case. The problem with this 
approach is that continuum value of J'b"'^ in fig.(|l|) is the extrapolation of the 
ratio of two quantities (in this case Zl and (fi^)'^/^, see eq.(|l^)) which have very 
similar functional dependencies on a. Any slight systematic effect in either the 

* as opposed to asymptotic scaling - see sec. 4 
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numerator or denominator will swing the ratio, which therefore could significantly 
affect the extrapolation. 

In this paper we show that such a systematic effect is present in the form of 
scaling violations for ^ 2 GeV. In the approach presented here, the functional 
dependence of the numerator and denominator are determined separately, and 
then compared to check for scaling. We then present a method which does not 
suffer from the above problem and enables a continuum estimate of /^*"* to be 
made. 

We begin by assuming a naive scaling of the quantity appropriate for 
ie. Zl (seeeq.dg)) 

Zl ~ e~^^fB/9\ (13) 

We will justify this choice of functional behaviour in the next section and show 
that, for our purposes, there is no loss of generality in eq.(^3]). Assuming this 
relationship, we now plot in fig.(@) log{ZL) against 1/ g^. Later we will fit this plot 
to a straight line to extract Sf^ Q. We notice immediately from fig.(0) that the 
data fall in a roughly linear band, and with a relatively small spread. One could 
imagine using this plot to check future calculations of at smaller a values. It 
is obvious that the behaviour of Z^ with g"^ apparent in fig.(0) is more clearly 
manifest than the behaviour of /^*"* with a (see fig.(|l])). 

We now determine S*/^ from the slope of log{ZL) against 1/g'^. For the 
Wilson data, we choose two intervals: 5.7 < /5 < 6.0 and 6.0 < [3 < 6.3. The 
motivation for these two intervals is that there is evidently a dependence of Sj^ 
on g'^. This will be discussed in detail in the next section. We obtain the values 
reported in the third column of table ^ For the clover data, we choose the interval 
6.0 < /? < 6.2 since the /f data is available only at /? = 6.0 & 6.2. 

Of course, the values of Z^ and their scaling with g'^ are not enough to 
determine the scaling of /^*"* itself. One needs to also study the g^ behaviour 
of some Jlf in order to set the scale. Four such choices are y/a, Mp, and the 
IP — IS splitting in charmonium. In tables ^ & ^ a sample of the published 
data on from various collaborations is presented with the references cited in 
column 1. Note that for aT^, we have used the renormalization constants shown 
in the tables. Specifically these were obtained with the "boosted" perturbation 
theoryQ in the Wilson case [^, and the non-perturbatively obtained values in 



the clover case [27, 28 



^Note that for ease of presentation, we have pfotted values obtained with both Wilson 
and clover actions together on the same plot. A priori, we cannot expect the same value 
of Sfj^ from both actions, so in the determinations of S'y^ we fit results from the two actions 
separately. 

® In this case we choose the "boosted" coupling gy — /u^ where Uq is the average plaquette, 
see Eyl. 
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Ref. 



Action P = 6/g^ Zl 



ELC 


Wilson 


6.0 


0.22(2) 


2.0(2)'' 


0.89 


310(25)(50) 


WPC 
WPC 
WPC 


Wilson 
Wilson 
Wilson 


5.74 

6.0 

6.26 


0.543(20) 
0.231(10) 
0.125(8) 


1.119(8)*= 

1.88(2)*= 

2.78(2)*= 






FNAL 
FNAL 
FNAL 

"PM A T 


Wilson 
Wilson 
Wilson 
vv iison 


5.7 
5.9 
6.1 
u.o 


0.564(28) 
0.250(14) 
0.135(13) 
u.uyy (^o ) 


1.15(8)^ 
1.78(9)^ 
2.43(15)^ 


0.63^^ 
0.65^ 
0.68^ 

U.Oo 


271(13)(20) 
241(13)(13) 
215(21)(14) 
ooK(^ 7\('\ A\ 

Z,ZoyL 1 ) ) 


BLS 


Wilson 


6.3 


0.094(6) 


3.21(9)(17)'^ 


0.7^^ 


235(20)(21) 


UKQCD 
UKQCD 


clover 
clover 


6.0 
6.2 


0.21lt? 
0.1171^ 


2.0+f 
2.7tp 


0.78^ 
0.79^ 


^00-10-42 
oc;o+16+105 


Ken 


Wilson 


6.0 


0.184(7)"^ 


2.0 


0.70^ 


224t? 


APE 
APE 


Wilson 
clover 


6.0 
6.0 


0.23(3) 
0.23(3) 


2.11(5)(10)^ 
2.05(6)^ 


0.89 
0.89f 


350(40) (30) 
370(40) 


APE 


clover 


6.2 


0.111(6) 


3.0(3)^* 


0.81^^ 


290(15)(45) 



Table 1: Values for Zl and fg'^* from various group's work. 

°- this Zl value was obtained from eq.(|lT]), ie. the Z^ value was not exphcitly 

published 

^ from averaging scale obtained from and Mp 

from the IP-IS value at /3 = 6.1 using 1-loop asymptotic freedom to ex- 
trapolate to /5 = 6.3 

from /^r 

^ from the string tension, a 

^ oT^ from IP — IS* charmonium splitting 

g i^Ren fj^Qj-^ standard perturbative result 

h i^Ren fj^Qj^ boostcd, tadpolc improved analysis 

I ^Ren fpQj^ boosted analysis 
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Action Fitting interval Sj^ Sf^ Sip-is 



Wilson 5.7 < /? < 6.0 2.39(10) 2.06(3) 1.53(7) 1.8(2) 2.2(4) 
6.0 < /? < 6.3 1.57(11) 1.52(6) 1.35(13) 1.7(2) — 



Clover 6.0</3< 6.2 2.1(2) 1.52(6) 1.7(2) 2.1(3) — 

Table 2: Values for Si obtained from fits to the data in tables |^ and to 
egs.(|ig,[7^. Note l3 = 6/g^. 



Again we assume the data to have the following functional form (see eg. [^): 

= ^ ~ (14) 

where i = a, Mp, IP — IS", and we have used eq.(|). As in eq. ([T3|) , this func- 
tional form can be assumed with no loss of generality. 

In figs.(0-|^), log{ai) is plotted against 1/ g"^ for the Wilson data from table ^ 
From the gradient of this plot fitted in the same intervals as the fit of loglZ^), 
we obtain the values of Si, i = {a, Mp, f^,, IP — IS} reported in table ||. 

We now turn to a discussion of the continuum limit of /^*°*. The functional 
behaviour of /^*"*, using eqs.( [Tl| . [l3|jI^ ) is: 



fr{9')-^,=e-'^'^B-'^^/'\ (15) 

We have ignored the g"^ dependence of the renormalization constant Zj^"', and 
will justify this below. Thus, the issue of scaling is addressed in this analysis by 
a comparison of Sf^ and Si. We beheve this is a cleaner method of studying the 
scaling of /^*'^* since the systematics present in both Zl and can be isolated 
and studied. 

From table |^ for (3 = Q/g"^ < 6.0 there is a clear statistical evidence for a 
violation of scaling in the Wilson data, ie. the Si are not all compatible. 

For 6.0 < f3 < 6.3, in the Wilson case, all the Si are within around one 
standard deviation of the value 1.5. In fact both the S^j and Sf^ values in table 
1^ are probably a little too high because there is no or ajj data in this range 
above P = 6.2. The conclusion is therefore that /b"^*, a, Mp, and IP — IS* 
mutually scale within present statistical errors for ^ 2 GeV. This is the main 
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Ref. 

w 

i- 
m 

\m 
\m 

\§- 

m 

\m 

m 

m 

i- 

in 

\m 

m 





- WPC 

- FNAL 
GFll 

- WPC 

- WPC 

- WPC 

- FNAL 
GFll 

- WPC 

- APE 

- APE 

- FNAL 
GFll 

- WPC 

- APE 

- WPC 

- WPC 
ELC 



P = Q/g' 

5.7 

5.7 

5.7 

5.74 

5.8 

5.9 

5.9 

5.93 

6.0 

6.0 

6.0 

6.1 

6.17 

6.2 

6.2 

6.26 

6.4 

6.4 



1.025(3) 



1.272(6) 
1.55(2) 



1.88(2) 



^Ren 



2.55(1) 



3.38(11 



.42(2) 


1.25(5) 


0.75 


.44(3) 






.99(4) 


2.00(5) 


0.77 


.25(10) 






.23(5) 


2.21(8) 


0.78 


.18(9) 


1.97(8) 


0.78 


.77(4) 


2.82(7) 


0.79 


.88(24) 


2.96(24) 


0.79 


.69(32) 






.70(15) 


4.0(6) 


0.80 



1.15(8) 



1.78(9) 



2.43(15) 



Values for a ^ obtained from various group's work using the Wilson 



Table 3: 

action. 



Ref. 



[3 = 6/^2 ^ 



^Ren 

f-r: 



[|T|] - APE 
0] - APE 

m\ - APE 

UKQCD 
APE 



m 



m 



6.0 2.05(6) 2.11(11) 1.09(3) 

6.0 1.92(11) 1.94(5) 1.09(3) 

6.0 1.95(7) 1.78(9) 1.09(3) 

6.2 2.7(1) 3.2(2) 1.04(1) 

6.2 3.05(19) 2.73(17) 1.04(1) 



Table 4: 

action. 



Values for a obtained from various group's work using the Clover 
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result of this paper. Results which suggest a decrease of f§°'^ as a — > can 
now be understood to be due to the inclusion of /5 < 6.0 data which suffer from 
systematic effects due to the lack of scaling of the quantities involved. 

In the clover case, there exists only two collaboration's results, and at only two 
(3 values, and so an interpretation of the results, in this case, may be premature. 



To elaborate, if the UKQCD results [16| are not included in the determination of 



the Si, we obtain the following values: 

=2.1(2), 5m, = 2.1(3), % = 1.7(3), (16) 

to be compared with the values in the last row of table ^ Evidently the Sm and 
Sf^ values interchange when the UKQCD results are not included. This suggests 
that more results, particularly at larger /? are required to settle the Si values in 
the clover case. In any case, the results suggest that f^^^,Mp and /^r mutually 
scale also in the clover case (for (3 ^ 6.0). 



4 Discussion 

In the rest of the paper asymptotic scaling, lattice artefacts and the continuum 
value of f^""* will be discussed. 

Asymptotic scaling is where the dependence of lattice quantities fif is the 
same as that predicted by weak-coupling perturbation theory. This behaviour 
can be obtained by integrating the beta-function of the theory, 

;,(,^)=„^ = 2Ajij + 2ft^ + 0(A (17) 



to obtain, 
where 

/?o 1 are given by: 



a = K-' fpTig'), (18) 
fpT = {9')^ e'W^, (19) 



(llA^-2nf) ^ 34^. 10,, (N^-1) , , 

is the number of flavours, the number of colours and A is some constant of 
integration. The subscript "PT" refers quantities obtained from (second-order) 
weak-coupling perturbation theory. For the quenched theory^, rij = 0, so we have 
from eq.(n) 



^see eg. p2i for a discussion of asymptotic scaling in the unquenched theory 
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where we have ignored higher order terms (ie. set (3i = 0). We expect, for small 
enough (7^, that Si = SpT- From table |^ we see that this is not the case. In the 
following, we discuss the possible causes for this discrepancy. 



4.1 Effects of higher order terms 

The first and most obvious explanation for the inequality. Si 7^ SpT, is that 
for the values of g"^ in table ^, the effect of higher order terms (/3i,/32, •••) is 
significant. Note moreover that the downward trend of the 5*^ towards Spt as g^ 
decreases is at least consistent with the prediction of weak-coupling perturbation 
theory and the declining importance of the higher order corrections in this limit. 

However, to study this hypothesis quantitatively, we first note that from 
eq. ([Tl|) , the lattice prediction of /^*"* is given by: 

3/3i 

fstat ^ 2Ren^g2^ Zz.(/) (/)^ e''--/^\ (22) 

where we have assumed the two-loop formula (eq.(|T8|)) for a, and have included 
a g^ dependence in Z^*^". Zj'^"' has the following expansion in g"^: 

Zf;"(/) = l-e/ + 0(/), (23) 

where e ~ 0.2. In this formula we have taken only the matching between the lat- 
tice and continuum effective theories and have ignored the anomalous dimension 



33] (see also [0). It can easily be demonstrated that this extra factor does not 



affect the following discussion. 

Since /^*"* is a fixed physical number, we have, using eqs.(|r3|, |22| & p3D 

-3/3i 

Z^(g^) ^ (1 + eg^) (/)^ ~ e-'^'fB^a'ya' (24) 



where we've now allowed for a g dependence in 5/^. Solving eq. (|2^) for Sfg{g ), 
using the fact that e is small, we obtain: 

SfA9')-^Spr-^f + y'. (25) 

Thus for e ~ 0.2 we get 

Sf,{g^) ^ 1.20 - 0.07/ + 0.02/, (26) 



which is not compatible with the values in table |[ More significantly, eq.(|26|) 
does not explain the strong dependence of Sjg on the g"^ interval (see table |^) 
or even the fact that Si > Spt- Thus this analysis suggests that higher order 
effects cannot explain the lack of asymptotic scaling of /^*"* (ie. Sfg ^ Spt) in 
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Fitting interval 'S^ 5*^ from eq.(p7|) 



5.7 < (3 < 6.0 1.39(2) 
6.0</?<6.3 1.26(5) 



1.07 
1.09 



Table 5: Values for 5*^ obtained from (i) fits to the string tension data in table 
\^ to eg. ^m), but with g ^ gv (column 2), and, (ii) a theoretical evaluation using 
eq.l^ (column 3). 



table @. A similar analysis of the affects of the higher order terms ui Si, i = 
{a, Mp, IP — IS"} leads to the same conclusion. The only difference in these 
cases is that the definition of the Si in eq.(|lj) includes the Zj^^"' factor, whereas 
the definition of Sj^ in eq.(0) does not. 

The derivation of Sfg{g^) leading to eq.(^5|) proves the statements above 
which stated that the functional form chosen in eqs.(|T^,|T^) is quite general for 
our purposes. Any g"^ dependence in eg. Z"^*^" can be factored into the definition 
of Si. It is also appropriate to note here that the contribution of Z^'^"' to Sfg{g'^) 
is ~ e/9 ^ 0.03 which is much smaller than the typical statistical errors present 
in Sfg{g'^) in table H. This justifies ignoring the g^ dependence of Zf^'^ in eq.(|l5D. 

Recently, effects due to higher order terms in lattice perturbation theory 
have been addressed 
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In this work, it has been suggested that the coupling 
constant appearing in the lattice action, is a poor choice of expansion parameter 
and the use of a "boosted" coupling, gy, was advocated. A typical choice of gy is 
9v ~ 9^ I'^ti where u\ is the average plaquette. A straightforward fit of Oo- as in 
eq. (|T^) , but with g replaced by gy., leads to the values of S^ in table |, column 
2. 



A theoretical evaluation based on that leading to eq.(^), but with g replaced 



with gy gives the result: 



s 



Pi 



PT 



1.20 



9v 



O.Olgl 



(27) 



e does not appear in eq.(p7D since it is zero for the string tension (ie. Z^^" = 1). 
The values obtained from this formula are shown in the last column of table ^ 
Again, the theoretical predictions do not match the results from the data. Also 
it is significant that the trend of S^ with g'^ cannot be reproduced, even with a 
boosted coupling. 



We do not choose to fit the data for a 



-1 

Mp,U,lP-lS 



to the boosted asymptotic 
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scaling prediction since the larger errors in these cases make the conclusive inter- 
pretation of results difficult. 



4.2 Lattice Artefacts 

The above discussion on higher order terms is an entirely continuum issue - it does 
not include any effect which is purely lattice in origin. In the following we discuss 
artefacts of the lattice formulation. Recalling the discussion in the introduction, 
we list the "lattice parameters" that are involved in these calculations^ { rriq, L, 
Quenching, a{g'^) }, where niq is the (light-) quark mass and L is the physical 
extent of the lattice in fermi. We discuss the effects of each of these parameters 
in the following to try to determine the cause of the observed lack of asymptotic 
scaling. 



"mg" effects 

The lattice values of /^*°* in this study are all the values obtained after an 
extrapolation in the light quark mass, ruq, to zero. Thus there is no problem 
associated with the light-quark mass not being adjusted to its physical value. One 
may, however, worry about a systematic error due to the extrapolation in ^ 
This is unlikely to cause a problem since the dependence on of Zl and indeed 
Mp, f.„, and many other quantities is mild (see eg. table 3 in Presumably 
also, any systematic effect associated with nig does not depend greatly on g"^ over 
the range studied in this analysis. 



Finite Volume effects 

The effects of finite L on many physical observables has been extensively 
studied P, For example, within present statistics, for L ^ 1.5/m, /^*"* 



is not a function of L [|T^ . This bound is not entirely satisfied by all the data in 
this analysis. So to study this explicitly, we take the values from table |l| at 
P = 6.0 and plot them in fig.(0) as a function of the L used for each simulation. 
There appears if anything to be a decrease in Zl, and therefore correspondingly 
/^*"*, as L increases, and this is contrary to the observed behaviour ||l3l (see 



also |T5[)- For this reason, and because generally speaking, different g'^ values 
in table |] have their lattice sizes chosen such that L f» const, we do not believe 
that finite volume effects are to blame for Si not being equal to its asymptotic 

^ In the past there was concern that /^*"* may also be dependant on the "smearing" size 
of the interpolation operator used to extract the matrix element However, it is now clear 
that it is not the case [§ |, |l|| . 

^Normally simulations are performed with niq ^ 100 MeV >> m„^d. 
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freedom value, Spt- 



Quenching 

The numerical effects of quenching on lattice calculations are normally diffi- 
cult to uncover, but in the case of Si this is not the case. To see this we note that 
the quenched approximation is really an effective theory of QCD where the sea 
quark loops are totally neglected and the coupling, g"^, is adjusted to try to com- 
pensate for these missing interactions. In an ideal world, this adjustment would 
be perfect, and predictions from the correctly-adjusted quenched theory would 
match those of the unquenched theory. However, due to the complicated structure 
of the interactions, this is presumably not the case! This means that, for example, 
instead of using the experimental (ie. unquenched) value of = Mp = 770 
MeV to set , we should use the value of Mp in a world without sea quark loops 
but where other physical quantities (such as a etc.) were as close as possible to 
their experimental values in our unquenched world. That is we should use 



j^Quenched— World 

< = -^^L (28) 

p 

cf. eq.(P). The difference between this prescription and the usual one (where 
instead of T^^Q^enc/ied-iuorid have Mp = 770 MeV) is simply a constant pre- 
factor, independent of g"^, and therefore this cannot affect Sm„- 



"O (a) -Effects" 

The final lattice artefact to be discussed is that due to the finiteness of the 
lattice parameter a. These so-called "0(a)" effects have long been studied (see eg. 



25, 36]) and are known to play an important role in matrix element calculations 
Z^, P^. They arise because the standard Wilson action replaces continuum 
derivatives by finite differences over nearest neighbours and is thus equal to the 
continuum action only up to 0{a). The clover action is improved to the extent 
that its predictions are correct to 0{a/log{a)) [0. (Note that the static theory 
on the lattice is also correct to this order |^.) The pure gauge sector of both 
lattice actions is correct to O(a^), and therefore the string tension, which is a 
pure gauge quantity, is correct to this order. 

Summarizing this discussion, we have for the Wilson case: 

af = nr + 0(a), (29) 
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and for the clover case: 



fif = fif^ + 0(a//o^7(a)), (30) 

where i = fs, Mp, IP — IS", and for the pure gauge case: 

fi^ = fif^ + 0(a^). (31) 

The fif ^ are the results that would be obtained with a perfect lattice action, ie. 
one correct to all orders in a. Using eqs.(^ & 0) together with eqs.(|^ & |T8|), we 
obtain, in the Wilson case for i = Mp, /^r, IP — IS, 



and, for i = a. 



a-\9') = K fpW) (1 - X^if^r), (33) 

Jpt['-) 

where Xi is the relative strength of the 0{a) (or 0(a^)) correction at g"^ = 1, 
ie. P = Q/g^ = 6.0. Here we have replaced api with fpT (see eq.(|T8D) since 
apL = A~^fpT up to exponentially suppressed non-perturbative pieces. 

Thus we see straightaway that in the limit where the effects of quenching are 
unimportant (and therefore Aj = Aj) a^^(g^) — aj^{g'^) = const, for fermionic 
Wilson quantities. In the future, with better statistics, this can be checked. 



We are now in a position to fit the data in table ^ to the appropriate eqs.(^ 
fc |33|). In this fit there are two free parameters: Xi and the coefficient, Aj. 

Taking all the a~^{g^) values in table |^ and fitting them to eq. (|33|) we obtain 
Xcr = 0.197(2) with a /dof of 11/4 which is shown in table ^ As a check we 
also fit a^^i^g"^) to eq.(p2D obtaining, as expected, a poorer jdof of 26/4 with 
Xg- = 0.343(4). This suggests that we are in fact isolating the O(a^) corrections 
in a/o", and that furthermore, the largish value of /dof in the fit to eq.(p3D may 
signal the statistical presence of even higher order terms (ie. of 0(a"),n > 2). 
The fit of a~^{g'^) to eq.(p3D is shown in fig.(§) as a solid line. In this figure 
we have also shown, as a dashed curve, a fit to the 2-loop asymptotic scaling 
formulae with the boosted coupling gy as discussed in sec. 4.1 (ie. eq.(|T8|) with 
g"^ ^ gy = g"^ /uq, where Mq is the average plaquette). As can be seen from fig.(^, 
the quality of this boosted asymptotic scaling fit is poor; in fact the /dof for 
this fit is around 10'^. On the other hand, the fit to eq.(^3D (ie. the asymptotic 
scaling formulae with an 0{o?) term) is very acceptable. This suggests strongly 
that the observed lack of asymptotic scaling in the string tension data can be 
simply explained by 0(a^ )-effects. 
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i 


a 






IP-IS 


fit to eq. 


m 


m 


m 


m 


X, 
Ai (MeV) 
Xl/dof 


0.197(2) 
1.780(5) 
11/4 


0.21(2) 
2.12(6) 

8/8 


0.31(5) 
2.4(1) 

8/5 


0.35(11) 
2.5(3) 
0.3/1 



Table 6: Values for Xi,Ai and x^/dof obtained from fits of the Wilson data in 
table to eqs. (^)^ as indicated. 



In the case of aMp,/^,ip-is Wilson data, a fit to eq.(|32D gives the values 

for Aj and Xj shown in table ^ These fits are shown in figs.(|9|-[llD. Again the 
/dof are very acceptable. 

Due to the quality of the fits we conclude that the most satisfactory explana- 
tion of the observed lack of asymptotic scaling (ie. S„^Mpj^,ip-is Spt) is 0(a) 
effects. This is the only explanation out of those discussed here which seems 
consistent with the data. 

The results for the A values in table ^ indicate that Act is significantly lower 
than the Wilson values for Ajv/pj^.ip-is, and that therefore results using will 
differ statistically in the continuum limit from those using a^^^ j-^ ip-is- However, 
a itself is in fact a poorly determined quantity since it relies on model calculations, 
and furthermore, from above arguments, the low value of A^- may signal the effect 
of quenching. 

For the clover case, we do not attempt to extract the coefficient of the a /log (a) 
term due to the fact that there are data at only two values of g"^ available. We 
await further data before attempting this analysis. Also, we choose not to perform 
a combined least squares fit of ^Mpj^^ip-is to eq.(^) together with to eq. (|33|) 
with a single A, since quenching implies that a single A is inappropriate. 

We can continue the analysis of the 0(a) effects by fitting as follows: 



^2/3--. ; ^^-Vpt(1)^- ^^^^ 

We obtain the values of X^^ and Xz^ in table |^. (Again we have fitted only the 
Wilson data.) Note that this functional form was chosen to mirror the fits of 
'^CT,Mp,/^,ip-i5 eqs.(|32| & ^). However, of course, in this case the physical value 
corresponding to Zj, (ie. fs) is not known, and there is nothing to set the scale; 
therefore A^^ is dimensionless. In fig.(|^) we plot the data against P (for the 
Wilson action) together with the fit to eq.(0) shown as a solid line. Note from 
this plot and from the relatively poor in table ^ that there appears to be some 
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fit to eq. p4| 




0.42(3) 




3.7(1) X 10-3 


Xl/dof 


39/9 



Table 7: Values for X^^, Xz^ and x^/dof obtained from a fit of the Zl values 
in table |7| to eq.j^S^). 



systematic effects remaining in the data. This is not entirely surprising due to 
the known difficulty of extracting Zl. 

We are now in a position to obtain a continuum value of /^*°*- Using 
eqs.(|n],p2|,p3|,pl), in the limit a 



/r (a = 0) = ^fr(^' = 0) (35) 

where in this limit, Zf^'^ig^ = 0) = 1 - '^§^ = 0.82 We use the average of the 
,/^,ip-is from table ^ ie. A = 2.2(2) MeV (where the error is statistical plus 
systematic combined in quadrature) and obtain, 

/|*'**(a = 0) = 230(35)MeV. (36) 

We take this as our best estimate of the continuum value of /^*°* from the lattice. 
Note that this is roughly equivalent to the values obtained from simulations at 
finite a values for ^ 2 GeV (see fig.(|l]), note though that this figure uses the 
scale from the string tension). Had we instead used A^. = 1.780(5), we would 
obtain /^*"*(a = 0) = 170(14) MeV. We do not prefer to choose this value since 
it appears that the values are contaminated by either quenched effects, or 
model dependences. 

Since the renormalization constant Zj^" between the lattice and continuum 
effective theories is evaluated at = 0, it is exactly determined. This is in 
contrast with the uncertainties in Zj'^"' at finite g^ which plague other approaches 
due to the uncalculated terms of order g^ and higher. 



5 Conclusion 

In this paper we have presented a method to obtain the continuum value of 
j!.^at fj^Qj-Q lattice data. This approach separately isolates the systematic errors 
coming from the dimensionless lattice quantity corresponding to /^*°* (ie. Zl) 
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and the lattice quantity used to determine the scale. These systematic errors 
and the lack of asymptotic scaling can be parameterised in terms of finite lattice 
spacing effects. Assuming this explanation a value of /^*°* = 230(35) MeV in 
the continuum limit has been obtained. In this method there is no error in the 
renormalization constant connecting the lattice and continuum effective theories. 
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Figure 1: Plot of f^""* against the lattice spacing as determined from the string 
tension, a^^, from various groups as listed in the legend. /^*"* has been determined 
using eq.(lJ) with a^^ = a'^^. See text for the definition of Zj^^ used. 
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Figure 2: Plot of log{Zi) from various groups as listed in the legend, g is the 
hare lattice coupling. 
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Figure 3: Plot of log {a ^) from the string tension. The reference is as appears 
the legend. 
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Figure 4: Plot of log {a ^) from Mp for the Wilson action. The references are as 
appears in the legend. 
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Figure 5: Plot of log{a ^) from for the Wilson action. The references are as 
appears in the legend. 
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Figure 6: Plot of log{a ^) from the IP — IS splitting in charmonium for the 
Wilson action. The reference is as appears in the legend. 
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Figure 7: Plot of Zl from various groups at (5 = = 6.0 against the spatial 
dimension in lattice units, L/a. The references are as listed in the legend. 
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Figure 8: Plot of a^^ from the string tension against (3 = 6/g'^. The solid curve is 
the fit to eg. j^SBj) (ie. the 2-loop asymptotic scaling formulae with an O(a^) term). 
The dashed curve is the fit to eg. ^TBj) with g"^ replaced by gy (ie. the "boosted" 
2-loop asymptotic scaling formulae). The references are as appears in the legend. 
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Figure 9: Plot of from Mp against (3 = Q/g^. The solid curve is the fit to 
eg. (ie. the 2-loop asymptotic scaling formulae with an 0{a) term). The 
references are as appears in the legend. All data is from the Wilson action. 
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Figure 10: Plot of a^^ from f^^ against (3 = Q/g'^. The solid curve is the fit to 
eg. (ie. the 2-loop asymptotic scaling formulae with an 0{a) term). The 
references are as appears in the legend. All data is from the Wilson action. 
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Figure 11: Plot ofa^^ from the IP— IS splitting in charmonium against (3 = Q/g'^. 
The solid curve is the fit to eg. (p^j (ie. the 2-loop asymptotic scaling formulae 
with an 0{a) term). The references are as appears in the legend. All data is from 
the Wilson action. 
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Figure 12: Plot of Zl against (3 = Q/g'^. The solid curve is the fit to eg. (ie 
the 2-loop asymptotic scaling formulae with an 0{a) term). The references an 
as appears in the legend. All data is from the Wilson action. 
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